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I. INTRODUCTION 


The radiation from slots and slot arrays in smooth, curved surfaces 
is pertinent to the design of flush-mounted antennas for high-speed 
aircraft and spacecraft. In this report the geometrical theory of 
diffraction (referred to simply as the GTD henceforth) is extended to 
treat the radiation from apertures in convex, perfectly-conducting 
surfaces. The GTD was introduced by Keller[l] as a systematic extension 
of classical geometrical optics to describe the diffraction phenomenon 
in terms of rays. Diffracted rays, similar to the ordinary rays of 
geometrical optics, are introduced to describe the propagation of the 
diffracted field; these diffracted rays are determined by an extension 
of Fermat's principle. Within the shadow region the geometrical optics 
fields are zero, and the diffracted rays alone account for the field. 
Diffracted rays are produced whenever an incident ray strikes a smooth 
surface at grazing incidence, or strikes an edge, corner, or a tip of 
a scattering object. At sufficiently high frequencies the excitation of 
the diffracted field is a local phenomenon; the diffracted field is related, 
to the incident field by one or more diffraction coefficients, which are 
functions of the boundary conditions and the local geometry at the 
point of diffraction. Conservation of power in a strip of rays on the 
surface is employed to deduce the behavior of the surface ray field; this 
introduces attenuation constants as parameters. Expressions for the dif- 
fraction coefficients and attenuation constants are found from asymptotic 
high frequency solutions to pertinent canonical problems or by boundary 
layer techniques. Away from the diffracting surface the behavior of the 
diffracted field along its ray path is the same as that of the geometrical 
optics field. The GTD is a high frequency approximation, which is asymptotic 
in the sense that the accuracy of the approximation generally improves as 
the frequency increases. 

Surface diffracted rays are also produced by apertures or slots 
in smooth curved surfaces; however, the excitation of these surface 
diffracted rays cannot be described by the curved surface diffraction 
coefficients alone as we do for the problem of scattering by a smooth, 
curved surface. It is therefore necessary to extend the GTD technique 
to treat the excitation of surface diffracted rays by apertures; we 
have achieved this extension by introducing launching coefficients 
which relate the excitation of the surface diffracted rays to the aperture 
distribution. If the field radiated from the surface is desired, the 
ordinary curved surface diffraction coefficients are used to determine 
the field of the rays shed tangentially from the surface rays; this 
field accounts for the radiation field in the shadow zone of the 
slot (aperture) in a smooth curved surface. In the region adjacent to 
the shadow boundary (which divides the region of space surrounding the 
antenna geometry into the illuminated region and the shadow region), the 
component of the field directly radiated from the source is represented 
by Fock-type functions. In the illuminated region the incident radiation 
from the source (this does not include the diffracted field components 
which may be present if the curved surface is closed) is treated by 


1 



geometrical optics. Usually, these solutions blend smoothly so that a 
continuous total field is obtained in the entire region of space 
surrounding the curved surface. 

Our reason for using the GTD method stems from the significant 
advantages to be gained; namely 

a) it is simple to use, and yields accurate results; 

b) it provides some physical insight into the radiation 
and scattering mechanisms involved; 

c) it can be used to treat problems for which exact 
analytical solutions are not available; and 

d) it can be extended in a straightforward manner to 
treat the mutual coupling between apertures. 

The extension mentioned in d) is achieved by introducing an attachment 
coefficient (related to the launching coefficient) which converts the 
field of the surface diffracted rays to the surface current (the field 
associated with the surface diffracted rays is not the field on the 
surface). This extension will be the subject of a later report. Even 
though the GTD is valid for high frequencies, it is surprisingly 
accurate in a number of cases for lower frequencies; the latter is 
shown to be true for some of the examples presented in this report. 

Also we note that as a ray optical technique the GTD is consistent with 
the reciprocity principle for electromagnetic fields. 

Sensiper[2] and Wait[3] have given asymptotic high frequency 
solutions for the radiation from slots in conducting circular cylinders. 
Belkina and Weinstein[4] have obtained asymptotic high frequency results 
for the radiation from slots in conducting spheres. However, none of 
these earlier results[2,3,4] were cast into the GTD form, and hence their 
solutions could not be directly extended to treat the radiation from slots 
in arbitrary convex conducting surfaces. One could use Fock theory 
and reciprocity to find the radiation field of apertures[5]; however, 
this method cannot be extended to determine the mutual coupling between 
apertures. Furthermore, the results given in [5] are not as accurate 
as the ones presented here. Balanis and Peters[6] have analyzed the 
radiation from slots in conducting cylinders by using the GTD; their 
approach is different from ours and is based on an approximate combination 
of wedge diffraction theory and surface ray theory. However, their results 
do not yield a continuous pattern across the shadow boundary unless they 
include higher order multiple diffraction between edges. Even with this 
included, the pattern still has some error near the shadow boundary, 
although it is continuous. 

Section II describes our extension of the GTD to analyze the 
radiation from apertures in arbitrary convex conducting surfaces. 
Expressions for the parameters and special functions involved are given 
in Section III. Applications of the theory to treat the radiation from 
slots in elliptic cylinders, spheres and spheroids are presented in 
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Section IV, where the numerical results obtained are compared with 
calculations based on rigorous solutions and with measurements. 

II. THE GTD FORMULATION 

The following development, based on the GTD, is heuristic. 

Physical arguments are used; ray trajectories are deduced from a 
generalized Fermat's principle, the variation of the field amplitude 
along a ray path is found from power conservation in a tube (or strip) 
of rays, and the diffraction and excitation of waves are treated as 
local phenomena. The validity of the GTD has been established through 
its application to a wide variety of antenna and scattering problems, 
and its present application is tested on certain canonical problems, 
which also yield the unknown parameters of the GTD solution. 

Consider the radiation from a slot or aperture in a curved, 
perfectly-conducting surface where the external medium is free space. ‘ 
The electric field at a point Q' in the aperture is E(Q'); it is assumed 
that the tangential component of this field is known. The magnetic surface, 
current associated with the aperture electric field is 

(1) Ks(Q') = E(Q') X n' 

A 

in which n* is the outward directed unit normal vector to the aperture 
surface at Q' . The infinitesimal magnetic current moment 

(2) dp^(Q') = l<3(Q')da 

radiating in the presence of the perfectly-conducting surface, which 
covers Uie aperture as well, is the equivalent source of the electric 
field dE(F) at P in the region exterior to the surface; here da is an 
area element at Q' . It will be shown how the geometrical theory of 
diffraction may_be used to calculate the incremental field dE(P). The 
electric field E(P) is then found by integrating over the aperture A 
shown in Fig. 1. 



Fig. 1. The equivalent source in the aperture A. 
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According to geometrical optics the region exterior to a perfectly- 
conducting surface with a magnetic current moment at Q' is divided into 
an illuminated region and a shadow region by a plane tangent to the 
surface at Q'. This plane is referred to as the shadow boundary. A 
portion of these two regions adjacent to the shadow boundary is a 
transition region. The angular extent of the transition region is of 
the order (kpg)"V3 in radians; pg is the radius of curvature of the 
surface at Q'^in the direction of the field point. The three regions are 
depicted in Fig. 2. We will describe how the high-frequency far-zone 
field may be calculated in the three regions. Let us begin with the 
field point in the shadow region (see Fig. 5). 

A. Shadow Region 

I 

The equivalent magnetic current moment Kgda' at Q' excites waves 
• which propagate along ray paths to the field point P; according to the 
^ generalized Fermat's principle proposed by Keller, the distance between 
Q' and P is an extremum (generally a minimum) on a ray path. With the. 
field point at P 3 in the shadow region, see Fig. 2, it is evident that 
a part of the ray path which minimizes the distance between Q' and P 3 
must lie on the perfectly-conducting surface, where it is a geodesic. 

In the case of simple surfaces such as the spherical surface, the 
cylindrical surface, the conical surface and the plane surface the 
geodesics are known and they are easy to describe; otherwise they can 
be found from the differential equations for the geodesic paths, which 
is a formidable but straightforward exercise; see Reference [25]. 

Now let us consider the behavior of the field associated with the 
surface ray (which we will henceforth refer to as the surface ray field): 

j(i>n-kt) 

(3) a(t) = A(t)e ° 

in which A(t) is real, t is the distance along the surface ray, 410 is 
the phase at t= 0 , and a time dependence of is assumed and sup- 
pressed. As it propagates along a geodesic on the curved surface, the 
energy flux associated with this field is proportional to A^dn, where 
dn is the infinitesimal distance between adjacent surface rays. It is 
assumed that the energy flux between adjacent surface rays is conserved. 
It follows from the generalized Fermat's principle that the surface ray 
sheds rays tangentially as it propagates; thus energy is continuously 
lost from the surface ray field, and it is attenuated. This may be 
expressed as 

(4) (A^dn) = - 2a(A^dn), 

where a is the attenuation constant. The above equation is readily 
integrated between t^ and t to give the surface ray field 
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ILLUMINATED REGION 
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Magnetic current moment on a curved surface showing the different regions of radiation. 



(5) 


a(t) = 


a(t„) 


dn(tQ) 

dn(t) 


-[jk(t-tQ)+ 


rt . 


a(t' )dt' ] 


This expression must be modified when there is a caustic due to a 
point source on the surface at Q‘ , where t=0. For to small driUo) 
tod4)Q» where di|;o is the angle between adjacent surface rays, see h'g. 3. 



Moreover, a(t) must be independent of to‘» hence lim a(to)|t^ exists 
as to^ and we define it to be C. It follows then that 


( 6 ) 


a(t) = C 




-[jkt + 


rt 


a(t')df] 


The constant C is proportional to the strength of the source at Q'. To 
facilitate the discussion Qf the vector source d^ni '•t is desirable to 
introduce the unit vector n which is normal to the convex surface at 
Q and the unit vector t which is tangent^to the surface ray at Q, see 
Fig. 4. A similar pair of unit vectors n]^, may be defined at Q'. 

In addition, the unit vectors 6=-nxt and b'=-n'xt' are defined at Q and 
Q', respectively. For the sake of simplicity, the following development 
is restricted to surface rays with zero torsion. 

The source d^m may be resolved into components parallel and 
perpendicular to tne surface ray of interest at Q' . For surface rays 
with zero torsion, it is assumed that these two components induce surface 
ray fields which propagate independently of each other along the 
geodesic arc between Q' and Q. First let us consider the surface ray 
field excited by the component of the source 6 • dpT^i which is 
perpendicular to the surface ray (S' = B when the torsion vanishes). 

From canonical problems, it is found that G • dp excites an infinity 
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TOP VIEW 



Fig. 4. Diffraction from a source Q' on a convex surface. 


of surface ray modes, and that the electric field of each of these modes 
is perpendicular to the surface. It follows from Eq. (6) that this 
electric field , 


(7) 


n • dEp(Q) 


% 



-[jkt + 



» 


where ap 
boundary 


is the attenuation constant of the mode for the hard 
condition, where the normal derivative of the field vanishes 
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at the surface.* The superscript h on a quantity denotes that it is 
associated with the hard boundary condition. The point Q is at a 
distance t from O' measured along the surface ray. The constant CA is 
proportional to b • d^m. The constant of proportionality is C LD, where 
lb is referred to as a launching coefficient and C is a conveniently 
chosen constant which is independent of p, so 

(8 ) C' = C lJ(Q') b . dF^(Q'). 

j. L 

We can now find the electrij field of the p^^ mode diffracted 
from Q toward P3; this field is n directed. From Fig. 4 it is apparent 
that the convex surface is the locus of one of the caustics of the 
diffracted rays. The other caustic is at a distance of p from the 
surface, as shown in the figure. From the geometrical theory of dif- 
fraction 


" • w ' ^pJ?feT 

where s is the distance between Q and P3 and Kp is a constant pro- 
portional to n • dTp(Q). In this case the constant of proportionality 
is Keller's diffraction coefficient for the high-frequency diffraction 
by a curved surface; hence 

(10) Kp = D^(Q) n • dEp(Q). 

In the present development we are interested in the far-zone 
field where s >> p. Combining the results of Eqs. (7), (8), (9), and 

(10) and summing over the modes of th£ surface ray field, the electric 
field at P3 due to the component of dp,|^ perpendicular to the surface 
ray is 

(11) 


in which 


(12) 

where drj; 


n . dE(P3) = C 6 • dpm(Q') F(Q',Q) 




-[jkt+ 


a"(t')df] 

0 P 


J d(j; 00 . . 

= dn/p is the angle between adjacent rays at Q. 


♦This approximation of the hard boundary condition for electromagnetic 
waves is adequate for the present discussion. A more complete treat- 
ment is given in Section IIIA. 
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Next let us consider the surface ray field excited by the com- 
ponent of the source t' • dpm parallel to the surface ray of interest. 
This electric field is again perpendicular to the surface ray, but also 
tangent to the surface. In this case it is clear that the surface ray 
field does not represent the actual field at the boundary, since it 
must vanish there. The precise nature of this field in the boundary 
layer is unimportant, because the field is merely a convenient artifice 
to relate the field at Q to the source. The surface ray field does 
not appear in the final expression for the field at P^. 

Following the development used to treat the 6' -component of the 
source, one obtains 


(13) 


b • dE(P3) = C t' 


dp^(Q*) 6{Q',Q) 



A 

as the b component of the electric field at P^, where 



condition, where the field vanishes at the surface. 

One may combine Eqs. (11) and (13) to_obtain the far-zone field 
at P3 due to the magnetic current moment dp^ at Q' 


(15) 




b G(Q',Q)] 



The attenuation constants, launching coefficients and diffraction 
coefficients which appear in the expressions for F and G are found 
from the asymptotic solution of canonical problems. The constant 
C also may be found from these solutions. The preceding discussion 
does not take into account the effect of surface ray fields which have 
encircled a closed surface a times. The inclusion of these fields is 
straightforward; one only needs to sum a power series in £, where the 
phase jumps at any caustics which have been crossed in an encirclement 
are taken into account. Since this phase jump equals tt/ 2 each time 
the caustic is crossed, it is easily taken into account by introducing 
a factor of j. Usually the contribution from the multiply-encircling 
waves is insignificant, unless the closed ray path is a few wavelengths 
or less. Also, in the case of a closed smooth surface the.re is generally 
diffraction from two points toward the field point; these are shown as 
points and Q2 in Fig. 5. 
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If the field point occurs at a caustic of the diffracted rays, 
so that infinitely many diffracted rays pass through P 3 , the geometrical 
theory of diffraction solution fails. However, it may be supplemented 
by an appropriate integral representation of the field, where the 
geometrical theory of diffraction is used to determine the equivalent 
sources of the integral representation. This is described in Appendix III. 

B. Transition Region 

The series representations for F and G in Eqs. (12) and (14) are 
rapidly convergent when the field point is deep in the shadow region. 
Usually only the first few terms are required to achieve reasonable 
accuracy, when the radii of curvature of the surface are larger than 
a wavelength or so. However as P 3 moves toward the transition region 
more terms must be added to maintain accuracy. As the field point 
moves into the transition region, the point Q approaches close to Q' 
and it is no longer desirable to treat the excitation, propagation and 
diffraction of the different surface ray modes separately. As a result, 
the series representations for the functions F and G are replaced by 
integral representations, and these functions are found to be pro-, 
portional to the Fock functions used to describe "creeping waves" in 
the asymptotic theory of plane wave diffraction by curved surfaces. 

From the viewpoint of reciprocity this result is to be expected; more 
will be said about this later. 

Also as Q approaches Q' , one notes that 


The ray configuration for a closed surface is shown in Fig. 5. 

The field of the ray Q'P 2 is calculated from the integral representations 
of F and G; on the other hand, ray Q'QiP2 is seen to have a considerable 
excursion on the surface and so the field of this ray is better calculated 
using Eqs. (12) and (14). 

C. Illuminated Region 


If the field point is at Pi in the illuminated region, the 
radiation from the source at Q' propagates along a straight line to P]. 
According to geometrical optics, the electric field at Pi due to the 
magnetic current moment dp^ at Q' is 

-jks 

06) E(P,) xdpjQ')|- ■ 

A 

where s is a unit vector directed from Q' to P-j , as shown in Fig. 6 . 
This equation may be written in a form similar to that employed when 
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the field point is in the shadow and transition regions, see Eq. (15); 
we proceed as follows. 



Fig. 6. An incident ray in the 
illuminated region. 


First let us define a plane of incidence which contains s and n', 
the unit vector normal to the surface at Q' . The unit vectors associated 
with the field, which are normal to the ray Q'P] and parallel and 
perpendicular to the plane of incidence are n and B, respectively. The 
unit vectors associated with the source, which are tangent to the surface 
at Q' and parallel and perpendicular to the plane of incidence are t' 
and 6, respectively. 

It follows from Eq. (16) that 

"jks 

(17) dE(P^) = C dp^(Q') • [b n 2 + t' b 2 cose] |- 

/V A 

in which cose = n'*s. From the equation above we see that 
(18a, b) F = 2 and G = 2cose 

according to geometrical optics. This approximation is found to be 
reasonably accurate deep in the illuminated region where 

e < f-- [kpg(Q')]‘'^l 


It joins smoothly with representations for the transition region, which 
are described in the next chapter. 
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In the case of the closed surface shown in Fig. 5, the fields of 
the rays emanating from Q' and subsequently diffracted from the points 
Ql and Q 2 are also included in the calculation of the field at P], These 
fields would be calculated from Eqs. (12), (14) and (15). 

It was noted earlier that ray optical approximations satisfy 
the reciprocity £rinciple. This can be confirmed for the solution just 
described. Let H(Q',P) be the magnetic field at Q' on the surface due 
to a unit magnetic current moment u at the exterior point P. This 
field can be found from an asymptotic solution of the diffraction 
problem posed by the perfectly-conducting curved surface. Let H(P,Q' ) 
be the m^netic field at P due to the differential magnetic current 
moment dp^^ at Q'. According to the reciprocity principle, 

(19) u . H(P,Q') = H(Q',P) • dp^ . 

It is found that expressions for the u-component of the magnetic field 
determined from the above equation are the same as those found from our 
application of the geometrical theory of diffraction. This checks the 
present development and shows that it is consistent with the reciprocity 
principle. 

D. Two-Dimensional Geometries 


The preceding development is easily modified so that it can be 
applied to 2-dimensional geometries involving infinitely long axial 
slots in the surfaces of perfectly-conducting cylinders. The tangential 
electric field in these slots may be axial or transverse Et, or a 
linear combination of the two. Thus two types of infinitesimal elementary 
sources at Q' are of interest: a line of axial magnetic current 

(20) dM^(Q') = + t' . dr(Q')df 

which is an ordinary magnetic line current flowing parallel to the z 
axis, or a linear distribution in the z-di recti on of magnetic current 
elements which are transversely directed, 

(21) dM^(Q') = b . dE(Q')dt' 


Note that b =+z here and dt' is an infinitesimal circumferential distance 
at Q'. In each case the upper sign is used if t is directed in the 
clockwise sense, the lower sign is used, if in the counterclockwise 
sense. 
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Many 2-dimensional problems of this type have been treated 
rigorously, so they may be used to check the accuracy of our GTD 
solution. Moreover, where asymptotic solutions exist, they may be 
used as canonical problems from which the parameters of the GTD 
solution can be found. Also such geometries are of interest, because 
their patterns are the same as those of finite axial and circumfer- 
ential slots in cylinders, when the patterns are calculated in the 
symmetry plane perpendicular to the axis of the cylinder. 

The surface rays emanating from the line source excitation are 
parallel so 


( 22 ) 



in Eqs. (12) and (14), and the amplitude variation along the ray shed 
tangentially from the surface ray to the field point is s-V2. in 
general 


(23) dM = zdM^ + t'dM^, 


consequently, the two-dimensional form of Eq. (15) is 

-jks 

(24) dr(P 3 ) = C dR.[bSF(q',Q) + t'b G(Q',Q)] S_ 

^ VT , 

in which the normalizing constant C may be expected to have a value 
different from that in Eq. (15). Bottj values of C are found from the 
solutions of canonical problems, and b = t z as noted earlier. 

III. THE PARAMETERS AND TRANSITION FUNCTIONS 

Expressions for the diffraction coefficients, attenuation con- 
stants and launching coefficients used in Eqs. (12) and (14) are 
given in this section along with the integral representations for 
the functions F and G employed in the transition region. The manner 
in which these integral representations transform to the geometrical 
optics field and the residue series used in the illuminated and 
shadowed regions is discussed. 

A. Diffraction Coefficients and 
Attenuation Constants 


The diffraction coefficients and attenuation constants depend 
on the local geometry of the surface, the wave number k, and the nature 
of the surface as described by the boundary conditions. Keller and 
Levy[7] have given the first order terms in the expressions for the 
diffraction coefficients and attenuation constants. However, before 
we present their results, it is desirable to examine the terms "soft" 
and "hard" boundary conditions mentioned in the preceding chapter. 
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This terminology is borrowed from acoustics. A soft boundary 
is one where the pressure field vanishes at the surface; it is also 
referred to as a Dirichlet boundary. On the other hand, a hard 
boundary is one where the normal derivative of the pressure field 
vanishes at the surface; this is also referred to as a Neumann boundary. 
Two types of surface ray modes have been assumed. For one type the 
electric field is in the binormal direction so that Ep = b Ep, and 
for the other, the magnetic field is in the bi normal direction sq _ 
that Hp = 6 Hp, and there is a normally-directed electric field n • Ep. 
The bi normally directed electric field clearly satisfies a soft or 
Dirichlet boundary condition at the surface, whereas the binormally- 
directed magnetic field satisfies the boundary condition 


(25) 



in which hb is the metrical coefficient associated with the unit vector 
b. The above boundary condition describes what we will refer to as 
a hard EM boundary. At high frequencies the second term is relatively 
small, so that the surface ray magnetic field satisfies a hard or 
Neumann boundary condition to a first approximation. Also Eq. (25) 
reduces to the hard boundary condition in the case of cylindrical 
surfaces where hb = 1. These observations concerning the boundary 
conditions are of importance in the paragraphs to follow. 


Let pg be the radius of curvature of the surface along which 
the surface ray is propagating in the plane containing the normal to 
the surface and the tangent to the surface ray. As mentioned earlier, 
Keller and Levy [7] have used first order asymptotic solutions for the 
diffraction of acoustic (scalar) and electromagnetic waves to deduce 
the attenuation constants and diffraction coefficients. For these 
canonical problems p = a, a constant. According to Keller and Levy 


(26) 


s _ 1 ,kaJ/3 
“op - I % 


(27) 


s 2 


for the soft surface, and 


(28) 


h 


a 


op 


= i(^)V3q 
a ^2 ' ^p 
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( 29 ) 


[D^ f 

op-^ 


'12 


. 1/2 


[Al(-qp)]^ 


for the 

hard surface. 

(30) 

AK-x) =1 

(31) 

Ai(-qp) = 0 

(32) 

Ai‘(-qp) = 0 

and the prime denotes 
the function. 


Voltmer[9] employing the same canonical problems as Keller and 
Levy, obtained attenuation constants and diffraction coefficients of 
improved accuracy by retaining higher order terms in the asymptotic 
solutions. Voltmer's results are presented in Table I, where it is 
seen that the corrections to the attenuation constants and diffraction 
coefficients are of order (2/ka)2/3. 

The first order approximations given by Eqs. (26) to (29) do not 
depend on whether the surface is cylindrical or spherical or on whether 
the wave is acoustic or electromagnetic. From Table I it is evident that 
this is no longer the case with the more accurate formulas. However 
an explanation of these second order differences is best accomplished 
after examining the high frequency diffraction by a more general 
surface, i.e., a surface of variable curvature along the ray path or 
of arbitrary curvature transverse to the ray path. 

Keller and Levy[15], Franz and Klante[21], Hong[10] and Voltmer[9] 
have considered the high frequency diffraction by general convex 
surfaces. Hong has obtained asymptotic solutions to the integral 
equations for the plane wave diffraction by a hard acoustic surface 
and the plane wave diffraction by a hard EM boundary. Voltmer has 
extended this work to soft boundaries, which are the same for acoustic 
and EM waves as we have noted. The solutions were carried out to 
second order, and they are functions not only of pg, the radius of 
curvature of the surface with respect to arc length along the ray 
trajectory, but also pg, pg and ptn> where the dot denotes a derivative 
with respect to arc length along the ray trajectory and p^p is the radius 
of curvature of the surface in the direction of the bi normal to the 
ray. Expressions for the attenuation constants are evident from the 
solutions; these are tabulated in columns C and D of Table II. On the 
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other hand, complete expressions for the diffraction coefficients can 
not be obtained from these solutions, because pg is assumed to be zero 
at the point of incidence on the surface, where the diffraction coef- 
ficient is evaluated. This condition was imposed to simplify the 
pertinent integral equations. The best available diffraction coefficients 
(more precisely, the diffraction coefficients squared) are given in 
Columns A and B of Table II. The incomplete portion of the second 
order term is indicated by (•••); it is a function of pg and pg.* In 
deriving these results it is assumed that pg/ptn > 1; furthermore as in 
the last chapter, it is assumed that the surface rays have no torsion. 

It is believed that the attenuation constants and diffraction 
coefficients listed in Table II are the best available at present and 
that they are adequate for most calculations, even though the expressions 
for the diffraction coefficient are not complete to second order. The 
improved attenuation constants are very important, because of the 
sensitivity of numerical calculations to errors in these parameters; 
corresponding errors in the diffraction coefficient are clearly less 
important to numerical accuracy. 

It is interesting to examine the general behavior of the 
attenuation constants and diffraction coefficients listed in Tables 
I and II. It is noted that all the corrections to the first order 
expressions of Keller and Levy are of order (2/kpg)2/3, and so their 
importance increases with a decrease in the. radius of curvature of 
the surface with respect to wavelength. All of the formulas in 
Table II reduce to those in Table I for the cylindrical and spherical 
surface, which is an important check on the asymptotic solution for 
the more general convex surface. The distinction between cylindrical 
and spherical surfaces in Table I is not required in Table II, because 
the term pg/Ap^p takes care of this difference; pg = p^p = a for the 
sphere, whereas ptn is infinite for circumferential rays on the cylinder. 
The formulas for the soft acoustic and soft EM surfaces are the same; 
this is expected because the boundary conditions are identical. On 
the other hand, the formulas for the attenuation constants and dif- 
fraction coefficients associated with the hard acoustic and hard EM 
surfaces are different, because the boundary conditions at the two 
surfaces are different. 

B. Launching Coefficients and Constants 

The launching coefficients have been defined in Eqs. (8) and (14). 

To determine the launching coefficients the radiation from a magnetic 
current nroment on a perfectly-conducting sphere and the radiation from 
magnetic current line sources on cylindrical surfaces have been 
analyzed. The asymptotic solution of these canonical problems and 


*The pg and J3g terms in the diffraction coefficient will be the 
subject of a^future investigation. 
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their ray optical interpretation is described in Appendix I, where it 
is found that for both cylindrical and spherical surfaces 


(33) H^^^‘{ka) D^, 

(34) = - j(jk|-)^/^ H^^^(ka) D^, 


where vp are the zeroes of the Hankel function in the first equation 
and the zeroes of the derivative of the Hankel function in the second 
equation. It is apparent that the relationship of the launching 
coefficient to the diffraction coefficient does not depend on the 
surface curvature transverse to the ray direction. For this reason, 
one may assume that 


(35) 2 


1 e 
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(36) 


-•IT 


LL' = 


e ^^(2TTk)^/2( 2 jl/ 3 ^.(_ ) 


P'L 


2 ^ 2/3 % 3 


2 tt ~i 




\t) 


J^P 


These expressions for the launching coefficients, where Dp and Dp are 
obtained from Table II, are the best available at present. It would 
be desirable to solve a canonical problem where the magnetic current 
moment or line source is on a surface of variable curvature as a 
further check. 


From Eqs. (A-26), 
that the constant 
( 


(37a, b) C = 


-jk 

4t7 


-ke 




j 


Birk 


(A- 36) and (A-59) in Appendix I it is found 
for 3-dimensional problems, see Eq. (15) 

for 2-dimensional problems, see Eq. (24) 


We note that C also could have been found from the geometrical optics 
expressions, i.e., from Eqs. (16), (17) and (18) for the 3-dimensional 
geometries, and from Eqs. (24) and (18) for the 2-dimensional geometries 
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C. Integral Representations for the Transition Region 


In the transition region which exists on either side of the 
shadow boundary, an integral representation for the field is necessary 
because the series representation for the F and G functions becomes 
poorly convergent; this integral representation yields a smooth trans- 
ition for the radiation field from the shadow region, where the surface 
ray modes are employed, to the illuminated region, where the geometrical 
optics field is used. As mentioned earlier, the angular extent of 
the transition region on each side of the shadow boundary is of order 
(kpg)"'/^. The expressions for the fields in the transition region 
are also deduced from the asymptotic solutions to appropriate canonical 
problems as described in Appendix II, and it is pointed out there that 
a first order approximation is usually adequate in the transition region. 

The integral representations for the field in the transition region 
can be expressed in terms of the well -tabulated Fock functions [22]: 

1 ) hard Fock function 


(38) 


g(^) = 


1 




^ Jr/2 


Wo(t) 


dT, 


2) soft Fock function 


(39) 


9 ( 5 ) = }= 


1 


e-JTC 

wj(Ty 


dt. 


which are described in Appendix II. 

Consider an infinitesimal portion of a surface ray on a perfectly- 
conducting surface of variable curvature. It follows from the. work 
of Fock[23], Wetzel [24] and Goodrich[5] that the variation of the field 
over this ray increment dt, and hence the variation of g(c) or g(^), 
behaves as though the surface were locally that of a cylinder or sphere, 
with a radius equal to that of the local radius of curvature in the 
direction of the ray p . Hence 


(40a, b) 


^ = 
dt 


g' (?) 


dt 


or 


If =9' (5) 


dt 


where 



(41) 


^ = 
dt 



9 


which is deduced from Eq. (A-71 ) after noting that adij; = dt. 

Generalizing F and G as given by Eqs. (A-79) and (A-92) in this 
manner for the 2-dimensional case 


(42) F = g(5) 




in which 
(44) 



follows from Eq. (41) and t is the distance between Q and Q', see 
Fig. 4. 


C, t > 0 when the field point is below the shadow boundary in the 
shadow region, and 5 , t < 0 when the field point is above the shadow 
boundary in the illuminated region, as shown in Fig. 7. When the 
field point is in the illuminated portion of the transition region, 
one visualizes the surface ray as traveling from its source at Q* to 
Ql , where it sheds tangentially back toward P. The Q' Qi P ray path 
does not obey the generalized Fermat's principle and therefore it is 
a pseudo ray system, but it does serve as a useful coordinate system to 
calculate the field at P in the illuminated part of the transition 
region. 


The function f(t) appearing in Eq. (43), is given by 


(45) 


f(t) =J 


1 


( 2 .1/3 


ft 


^ df 

0 Pg 2 


, -l<n'.s<0 

, 0<n'*s<l 


The introduction of f(t) into Eq. (43) gives us an algorithm for 
calculating the field in the transition region; furthermore this field 
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joins smoothly with the geometrical optics field deep in the illuminated 
region. The accuracy of the algorithm has been confirmed by applying 
it to a number of examples. f(t) reduces to Eq. (A-93) for surfaces of 
constant curvature, i.e., the cylinder and sphere. 



Fig. 7. An equivalent ray system for the field calculation 
in the illuminated part of the transition region. 


In the case of 3-dimensional geometries 
r 

^ -ikt 

^g(c)e 


(46) 


F=1 


diij 


g(s)e 


-jkt 


t ^ 0, 
t < 0, 


(47) 




■i ^1/3 

d^^^ ^kpTET^ 


G = S 




shadow region 
illuminated region 

t ^0, shadow region 
t < 0, illuminated region 


where ^ is given by Eq. (44), and the factor di|;o/dip is defined in Section 
I, see Figs. 3 and 4. In Fig. 4 recall that pd 4 ) = dn. This factor 
measures the change in angle between adjacent surface rays as they 
propagate from their source point to the point of diffraction Q] g. 

If dipo/dif; > 1 the surface rays converge and if d\i)o/d\p < 1, they diverge. 
di(;o/di|j is a function of both the ray and surface geometries; it may be 
found by using differential geometry. 
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At first glance it may seem strange thatxfdiji^/^ is not included 
in the expressions for the illuminated region. This may be explained 
with the aid of Fig. 7, where it is seen that the surface divergence of 
the rays emanating from the source to the point Q] is compensated by 
the subsequent convergence of the rays shed from Q, to the field point 
P in the far-zone. 

Equations. (42), (43), (46) and (47). also may be joined in the 
transition region with the corresponding surface ray mode expressions 
for F and G, see Eqs. (12) and (14). Noting that the former equations 
are a first order asymptotic approximation, whereas the surface ray mode 
expressions are approximated roughly to second order here, one would 
expect small differences in the values calculated from the two repre- 
sentations in their common domain of validity; however these differences 
should not exceed a few tenths of a dB in most cases. Deep in the 
shadow region the surface ray mode expressions are definitely to be 
preferred. 

The following summary of formulas is presented for the con- 
venience of the reader. It also facilitates the discussion in the 
next section. 

4 = [kpg(q')r’^^ 


in which Pg(Q') is the curvature of the surface at the source Q' in the 
direction of the surface ray. 


( 2 ) 

(15) 


(a) Three-dimensional problems ; 
d%(Q') = rE(Q') X n']da', 

- ^ 


dE(P) = ^ dpjQ') • [b n F + t' b 6] I 


,-jks 


(b) Two-dimensional problems : 


(23) dM = z dM^ + t' dM^, 

(24) dF(P) = dM • [b n F + t' b G] I: 


N^irk 


e 


-jks 


where b = 


+ z. 
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TABLE III 



F 

G 

1. 

Ft as in Eq. (18a) 

G] as in Eq. (18b) 

0 £ 0 < ^ - A 



2. 

?2 as in Eq. (46); 

G2 as in Eq. (47); 

^-A<0<^+A 

also Eq. (44) 

also Eqs. (44), (45) 

3. 

F^ as in Eq. (12); 

G^ as in Eq. (14); 

i+A<0'<'iT~<S 

also Table II and 

also Table II and 

c 

Eq. (36) 

Eo, (35) 


In application to two-dimensional problems 



6=0 unless an axial caustic of the diffracted rays is present in 
which case one employs Eq. (A-104) in the region TT - 6 < e ^TT» where 
6 = (ka)"l/3 in this inequality and a is defined in Fig. 5A. 
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IV. PATTERN CALCULATIONS 


In this Section, the expressions developed in the previous 
secions are employed to calculate the radiation from slots in cylinders, 
spheres, and spheroids. In dealing with the radiation from slots in 
perfectly-conducting cylinders, we formulate the problem in terms of 
its equivalent 2-dimensional representation; as noted earlier the 
radiation patterns of 2-dimensional, infinitely-long slots in cylinders 
are the same as those of the finite axial or circumferential slots 
when the patterns are calculated in the symmetry plane perpendicular 
to the axis of the cylinder. In the case of a slot in a perfectly- 
conducting spheroid, the slot dimensions are assumed to be infinitesimal; 
however this is not an essential limitation of the method, since an 
extended slot can be regarded as a continuous array of infinitesimal 
slots. Furthermore, these examples have been chosen so that the surface 
rays excited by the slot traverse geodesic paths with zero torsion;* 
this restriction is pointed out in the theoretical development of 
Section II. The radiation from slots in cylinders is treated first. 

The expressions for F and 6 may be determined with the aid of 
Table III on page 25. 

A. Radiation from Slots in Cylinders 

A perfectly-conducting convex cylinder excited by a magnetic line 
source on its surface Q' is depicted in Fig. 8a. Two kinds. of magnetic 
line sources Mz and Mt are considered; they are the equivalent sources 
of an infinitesimal axial slot, whose aperture electric field is in the 
circumferential and axial directions, respectively. Sources of this 
type have been discussed in Section ll, Eqs. (20) and (21), and in 
Appendix I, part 3. 

The OTD formulas for the far- zone field of this radiating system 
will be given in its shadow, transition and illuminated regions. The 
phase center of the field is chosen to be the origin 0 depicted in 
Fig. 8. 


*The modifications required of the present formulation when the surface 
ray possesses torsion (b';^b) is not well understood at this time; 
however Eq. (15} may be applied formally to torsional surface rays by 
replacing B by b' in its first term. In a number of examples it has 
been found that the resulting accuracy is acceptable for practical 
purposes; however this can not be asserted for the general case without 
further work. 
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Fig. 8b. 


(b) 

Ray geometry associated with a slot on a cylinder when 
P is in the transition region. 
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S2 


Fig. 8c. Ray geometry associated with a slot on a cylinder when 
P is in the illuminated region. 


1 . Deep shadow region 

The contributions to the far-zone field at P come from the 
diffracted rays shed tangentially from the surface at Qi and Q 2 toward 
P, with the result that 


( 48 ) 


-ke 


i- 

J/i 




jkr,*s jkr^*^ 

F3(Q',Q^)e ' +F3(Q',Q2)e ^ 


j 

JT , 
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i- 

J 4 

(49) ^ M 

^ 24 ^ ^ 


^ + G3(Q'.Q2)e'^''''' 


.•slp-jks 


where the position vectors of Q] and Q 2 are r-] and r 2 > respectively, 
s is a unit vector in the direction of P, and the approximations 


(50a) S.J s - r-| • s, 

(50b) S 2 s - ?2 * s. 


valid in the far zone have been employed. 


The contribution from rays which encircle the cylinder one or 
more times before sheding tangentially from Q-| or Q 2 can be neglected 
for all but the smallest cylinders with the hard boundary condition. 
There it can be taken into account by multiplying each term in the 
series of Eq. (12) by , 

r I u 


jknT-n 


(51) 


I e 

n=0 


a' (t')dt' 

Jo P 


■1 


-jkT- 


a"(f)dt' 


JoP' 


in which T is the circumference of the cylinder. The contributions from 
the field of these multiply-encircling rays were found to be negligible 
in the calculated patterns shown in this report. 

2a. Transition region (shadow part) 

Again the contributions to the field at P come from the diffracted 
rays shed tangentially from the surface at Q-j and Qp toward P, as shown 
in Fig. 8b, but in this case Qi is so close to Q' that the series of 
surface ray modes associated with this surface ray is poorly convergent 
and should be replaced by a Fock representation; hence 


i- 

(52) E = — 


_ A 

j k r • s 

F2(Q',Qi)e ^ +F3(Q',Q2)e 


jkr2*s1 e-jks 
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(53) 


-ke 




r jkr,-s 

Mt|_G2(Q'.Qi)e ' 


jkr„*s' 

+G3(Q'.Q2)e ^ 


^jks 

JT 


2b. Transition region (illuininated part) 

The contribution to the field at P comes from the ordinary 
diffracted ray shed tangentially from the surface at Q] toward P plus 
the contribution of the pseudo diffracted ray shed tangentially from 
Q2, as shown in Fig. 8c. The far-zone fields of the two line sources 
are given by Eqs. (52) and (53), but it should be noted that t in 62 
is the arc length Q'Q-i, which is measured negatively, thus 5 too is 
negative. 

3. Deep in the illuminated region 

Referring to Fig. 8c the dominant contribution to the electric 
field at P comes from the geometrical optics field, which propagates 
along the direct, straight line path between Q' and P. For the two 
magnetic current line sources being considered this is 


P.a p-Jks 


(54) e9-°* I 


JI7k 


JT 


(55) 


:9*0. 


J4 


J^k Js 


From the examples we have treated, it has been found that the 
geometrical optics field is so closely approximated by the field of 
the pseudo ray Q'Q-|P 


(56) 




F.(Q',Qi)e 

2j^ ^ ' 


jkri.s g-jks 

Js 


(57) 


5* _ “ke 


J4 


jkri-s g-jks 
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that in many instances the above equations can be used in place of 
Eqs. (54) and (55) throughout the illuminated region. 

In addition to the geometrical optics field there are much 
weaker contributions from the fields of the diffracted ray shed tan- 
gentially from Q 2 and the diffracted ray shed ^ngentially from Qi 
after traversing the surface in a clockwise direction from its source . 
Q'. These diffracted ray contributions are significant only for the 
small cylinders excited by the Mz line source {the hard boundary case). 
In summary deep in the illuminated region 


(58) 


E = E^*°*+ 


-ke 

2^2 


1- 

J 4 


+ F3(Q',q,) 

7T K L 


jkr 


•n^jks 


(59) = e9-°- . 


The patterns of axial slots in circular and elliptic cylinders 
have been calculated using the preceding formulas. In each case the 
electric field in the slot is circumferentially directed. Some patterns 
of an infinitesimally thin axial slot on circular cylinders of varying 
sizes are shown in Fig. 9. These patterns were calculated from Eqs. 

(48) and (52); it was unnecessary to use the geometrical optics repre- 
sentation deep in the illuminated region. In Fig. 9 the GTD patterns 
are compared with those calculated from the exact (eigenfunction) 
solution, and it is seen that for ka > 3 there are no significant 
differences between the patterns baseT on the GTD solution and the 
patterns calculated from the exact solution. Even for a cylinder as 
small as two wavelengths in circumference the error in the GTD pattern 
i s seen to be smal 1 . 

The pattern of a narrow axial slot in an elliptical cylinder is 
shown in Fig. 10. Eqs. (48), (52), and (56) were employed to calculate 
the GTD pattern; the finite width of the slot (0.34 wavelength) is taken 
into account by an array of 5 line sources in the aperture. This 
pattern is compared with a measured pattern of the same configuration 
in a plane perpendicular to the axis of the slot and cylinder[6]. The 
agreement between the measured and calculated values is good, since 
the discrepancy between the two halves of the measured pattern about 
its plane of symmetry is approximately the same as the differences 
between the measured and calculated values. This agreement is grati- 
fying because it confirms the accuracy of the GTD solution in an 
example where the surface curvature is variable. 
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RELATIVE AMPLITUDE 



<^ = AZIMUTH ANGLE (DEGREES) 


Fig. 9. Radiation patterns of a thin axial slot in a 
perfectly conducting circular cylinder. 


The GTD patterns of infinitesimal axial slots at positions of 
maximum and minimum radii curvature on an elliptic cylinder are shown 
in Fig. 11. As one would expect, the level of the pattern in the 
shadow region is much lower for case A, where the elliptic cylinder 
tends to block the rearward directed radiation more. 
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Fig. 10. Pattern of an axial slot on an 
elliptic cylinder. 


B. Radiation From Slots in Spheres and Spheroids 

In this section the pattern of an infinitesimal slot in a convex, 
finite surface of revolution is considered for the case where the slot 
is located on the axis of revolution. However, numerical results are 
presented for the spheroid and spherical surfaces only. 

The equiv^ent source of an infinitesimal slot is the magnetic 
current moment p^ located at the position of the slot and directed 
parallel to its axis. A densely packed array of such magnetic current 
moments can serve as the equivalent source of slot of finite size; 
alternatively it could be represented as a continuous distribution as 
in Eq. (2). 
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Now consider the radiation from an infinitesimal slot, with the 
magnetic current moment Pm = xp^, as its equivalent source, located 
on the axis of a convex body of revolution at z = Zg, as shown in Fig. 
12. From the symmetry of the problem, it follows that the radiation 
pattern can be completely specified in terms of the electric field 
pattern En(s,e) in the meridional plane parallel to the slot and the 
electric field pattern Ejs,e) in the meridional plane perpendicular 
to the slot; in other words, in terms of the H- and E-plane patterns. 





Fig. 12. Infinitesimal slot on the axis of a convex 
body of revolution. 


Introducing the spherical coordinate system (s,e,(f)) in the far zone 


(60) E(s, 0 ,(|)) = - $ E„(s, 0 )cos<{> - 0 Ej^(s, 0 )sin({). 

The GTD solution for the H- and E-plane patterns parallels the analysis 
carried out in the preceding section for the patterns of cylinders 
excited by the M. and M magnetic line sources, respectively, with some 
differences: ^ 


1 ) 

2 ) 


the surface divergence 



must be determined. 


the product of C and the radial dependence of the far- zone 
field is -jk e-jks/4Trs, 
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3) there may be a caustic on the ray path, 

4) there is a caustic of the diffracted rays at e = n for the 
configuration shown in Fig. 12. 

Figures 8(a,b,c) may be used to illustrate the pertinent ray trajectories 
for the body of revolution by considering the closed curve depicted in 
each figure to be the cross section of a body of revolution whose axis 
passes through the source point Q'. 

1 . Deep shadow region 


(61) 


r - ~ jk _ 

" 4ir "m 


G3(Q',Qi)e ^ +jG3(Q',Q2)e ^ 


■jks 


e^ 
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(62) Ej_ = 


- -jk 


47t ^m 


jkr,*s jkr^'S 

F3(Q',Q^)e ' +jF3(Q',Q2)e 


,-jks 


s 


except in the region of the axial caustic of the diffracted rays, 
TT-A<6<Tr, where Eqs. (A-105), (A-106) ar^ employed. 


2 . 

(63) 


Transition region 


E = D 
" 4tt ^m 


jkr, »s jkr«*s 

G2(Q',Qi)e ' +jG3(Q',Q2)e 


e-jks 

s 


has the same form as E„ except that F replaces G. Also note 
the comments in the preceding section concerning the differences in 
the ray paths and the calculation of F and G resulting from the field 
point being in the shadow or illuminated parts of the transition region. 

3. Deep in the illuminated region 


(64) 


E = -J— n 
‘"'I 4 tt 


2(z.s)ej''’"'*^ + jF3(Q',Q2)e 


jkr2*s 
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(65) 


•-L = Pm [ 


2 + jG,(Q',Q,)e 


jkr2*^ 


jkr, *s l_-jks 

+ JG3(q',Q,)e Jf 


where r' = z z . 

0 

In evaluating F and G the surface divergence factor must be 
determined. For the axially-excited bodies of revolution under study 
it is a simple exercise of differential geometry to find dij; /dt|;. 





Fig. 13. Diffracted ray geometry. 


Referring to Fig. 13 

dtij 

(66a) 

dll' Ptn^^"'*'! 
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where is the radius of curvature of^the surface in the direction 
perpendicular to the ray, ij^-] = cos"^(n*z), and p is the caustic distance 
at the Q-| , where the diffracted ray sheds tangentially toward P. Since 

(66b) p = p^j^ tani|;y, 

^ ■ "Eosli;^ * 

E |,and Ej^ given in Eqs. (61) through (65) may be substituted into 
Eq. (60) to calculate the far-zone field at any aspect. The preceding 
expression for E,, and Ej^ have been used to calculate the patterns of 
infinitesimal slots on perfectly-conducting spheres and spheroids. In 
Figs. 14 and 15 the patterns in meridional planes parallel and 
perpendicular to an infinitesimal slot are shown in Figs. 14 and 15, 
where they are compared with exact values calculated by King and Wu[12]. 
Although King and Wu actually calculated the currents on a perfectly- 
conducting sphere illuminated by a plane wave, their results can be 
compared directly with ours through the reciprocity principle. The 
GTD patterns are seen to agree very well with the rigorously calculated 
patterns; as a matter of fact in Fig. 15, they are seen to be in better 
agreement with these patterns, than a pattern calculated from a formal 
asymptotic solution given by Belkina[4]. 

The GTD patterns of axially-positioned, infinitesimal slots on 
oblate and prolate spheroids are shown in Figs. 16 through 19. No 
rigorously calculated or measured patterns were available for com- 
parison in these cases; however, the pattern for each spheroid is 
compared with the pattern of a slot on sphere with the same radius of 
curvature as the spheroid at its axis of revolution, where the slot 
is positioned. It is seen that spheres and spheroids with the same 
radii of curvature at their slot positions have essentially the same 
patterns in the illuminated and transition regions; the significant 
differences occur only in the deep shadow regions of their patterns. 
Considering the nature of the GTD solution described in Sections II 
and III, this result is not entirely unexpected. 
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GEOMETRICAL OPTICS 
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Fig. 14. Pattern of an infinitesimal slot on a 
sphere in a plane parallel to the slot. 
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Fig. 18. Pattern in plane n slot. 





V. CONCLUSIONS 


Keller's geometrical theory of diffraction has been extended so 
that the radiation from apertures (slots) in perfectly-conducting 
smooth curved surfaces can be treated. Configurations of this type 
are of interest in the design of flush-mounted antennas for aircraft 
and spacecraft. In extending the 6TD, a new parameter has been 
introduced, the launching coefficient, which relates the amplitudes 
and phases of the surface ray modes to the equivalent source in the 
aperture. This extension of the 6TD has been successfully supplemented 
in the transition region and at axial caustics so that the pattern 
can be calculated at all aspects. 

Although the present formulation is restricted to torsionless 
surface rays, it can be formally applied to surface rays with torsion, 
see Reference [25]. Examples considered in this reference suggest that 
moderate torsion does not result in excessive error in the pattern 
calculation. The propagation along surface rays with torsion is 
currently under study. 

Although the calculation of far-zone fields is described here, 
the method can be applied to calculate near-zone fields without any 
additional assumptions. Furthermore, by introducing a second new 
parameter, the attachment coefficient, the GTD can be extended to 
determine the tangential magnetic field at the perfectly-conducting 
surface, so that the mutual coupling between apertures can be calcu- 
lated; this will be the subject of a future report. 
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APPENDIX I 

DETERMINATION OF THE LAUNCHING COEFFICIENTS 


The canonical problems chosen to find the launching coefficients 
are the radiation from a magnetic current moment on the surface of a 
perfectly-conducting sphere and the radiation from two types of 
magnetic line sources on the surface of a perfectly-conducting cylinder. 
In the case of the first problem the asymptotic solutions for the E- and 
H-plane patterns are sufficient to determine the launching coefficients 
for the hard and soft boundary conditons. 

The radiation from a magnetic current moment on the surface of a 
perfectly-conducting sphere may be formulated conveniently in terms 
of a dyadic Green's function. The properties of dyadic Green's functions 
have been described by C.T. Tai in Reference [19]; the dyadic Green's 
functions for simple geometries, including the sphere are also given in 
this reference. 

Consider a perfectly-conducting sphere of radius a, whose center 
is at the origin of a spherical coordinate system as shown in Fig. lA. 


z 



Fig. lA. Infinitesimal slot on a sphere. 
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The electric field_of the magnetic current moment ^ p 
on its surface at R' = az, i.e., R' = a, 6 = 0 , is 


(A- 1 ) E(R) = p_ X • [V X G, (R',R)]|_ . 

' R'=az 

where 

( A-2 ) f ^ ( R‘ , R) = \ ( R' . R) + f 1 5 (R* , R) , 


(A- 3 ) G„(R' ,R) = - ( 2 - 6 „) ^ X 


m' (k) (k) + n' (k) (k) 

, oHin ^mn ^mn ^mn 


(A- 4 ) G,^(R'.R) ( 2 . 5 ^)^!^, 


a M^2)' ^(2) ^ ^(2)’ ^(2) 

I „nin „mn 

Loo 0 00 0 


(A- 5 a) , 9 p = - jn(ka)/h^^^(ka). 


(A- 5 b) bg = - [PaJn(Pa)] (pa^^ * 

o" 


(A- 5 c) 


Pa = ka. 


(A- 5 d) 


^0 


1 , m =0 

0 , m?^ 0 . 


posi tioned 


(k) , 
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(A-5e) R > R' = a. 


The time convention used by Tai is whereas our time convention is 

eJf^t; hence the superscript (1) in his formulas must be replaced by 

(2), which denotes a spherical Hankel function of the second kind, also 

i in his formulajs must b£ replaced by -j in ours. The spherical vector 

wave functions M and N are given in Chapter 11 of Reference [19], 

®mn ^mn 

0 0 

It is readily shown that 

(A-6a) X • Np (k) = X • rT , (k) = - in(n+l) , 

»mn 


and 

(A-6b) X • M (k) = X • = j n(n+l) jp(pg) . 

^mn 

0 


Substituting the above into Eq. (A-1 ) 
. 2 

(A 7l y 2n+1 

(A-7) E(R) - J^^T(ThT) 


f^PaJ’n^Pa^^ j'n^Pa^ I^Pa^n^^ (pa^^ K tt{2) 

"^eln 


'n'^a' 


"an '^a' 


L. \ 




r • * '' 

J.. , ^ ^Pa^n^Pa^^ u(2)/ ^lT^(2)/,.^ 

"i^n^Pa^ " r u(2)/ nh' \ ^PaM ^oln^"^^ 




by 


In the far-zone the vector wave functions may be approximated 


(A-8a) = 


m( 2) = .n e 


-jkR 
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in which 


p'(cose) . dpj(cose) 

(A-9a) s1n*e - gj cos+*. 

• _ dp'(cose) . p'(cose) 

(A-9b) RX = 35 sin*e •. ^335 cosM • 


Employing Eqs. (A- 8 ) and (A-9) in Eq. (A-7), the electric field 
in the far-zone is 


....2 


(A-10) E(R) = 


jk «-jkR ” - o_.i 

■n\ _ ^m e r /,-\n <:n+l 


me r e-ii-ri 

4tt . kR n(n+l ) 


where 


J^(ka) 

ka 


j„(ka) H^^^'(ka) 

3„(ka) j‘(ka) Hj[^'(ka) 

H<2)'(ka) 




f R X 


(A-lla) Z„(kr) =J^Z , (kr). 




(A-llb) = kr Q Z -j (kr). 


in which ® solution of Bessel's differential equation of 

order n+1/2 and Q is an operator defined by 


(A-12) 
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Substituting Eqs. (A-11) and (A-12) into Eq. (A-10) 


(A-13) E(R) ■v 


r /.xM £.n-ri 

n(n+l) 

n=i 


J ^(ka) 


n+7T 


,Q J i(ka)--rJ QH^^}(ka) 

H^"^j(ka) n+^ J 

n4^ 


'' dp|^(cOS 8 ) 


de 


COSdxf) + 




Q J i(ka) 




dP^(cose) 


V"-2 


Q H^^.j(ka) n-t^ 


sin(|)6 


pj(cose) 

where terms involving . ^ ^ have been neglected, because the solution 

Ka S 1 n 6 

for kasin 0 >>1 is of interest. 

The relationships 

(A-14a) ^ Pj(cose) = (-1)"" n{n+l) ^n^ 


and 

(A-14b) 


J 1 (kr) = ^ 


H^^^kr) + H^^j(kr) 
n+^ n+^ 


are employed in the development to follow.. 

1. H-plane (cf) = 0) Soft EM boundary 

Using Eq. (A-14), it follows from Eq. (A-13) that 
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(A-15) 


^ k p ^-jkR . 
T « 2 e 

t — (]> 


4 kR 


I (-l)"(2^)1q H<'j(ka) 
n=l ^ 1 n4^ 


H^^}(ka) 


nf^ 


Q H^^](ka)l^ 


dP”^(-cose) 


H'"{(ka) n+^ I 


The Watson transformation facilitates the asymptotic approximation 
of Eq. (A-15) and its transformation to a ray optical form suitable for 
calculating the field in the shadow region. Accessible descriptions of 
the method are given in References [16,20], Applying the Watson 
transformation to Eq. (A-15) and evaluating the resulting integral by 
the residues of its integrand at v = v , 


(ni6l r • 

(A-16) - 2 R I, -j2v • "f 

P ' 1 +e P 


H<l>(ka) 




i I 


(ka) 


. Q H^^^ka) 


,-l 


3Pvp_i/2(“Cose) 

39 


( 2 ) 

with vp being the complex orders for which '(ka) = 0, p = 1,2,3*** , 
For large ka P 


(A-17) Vp = ka + (^)^/^ qp e 

= ^4’ ^ ■’p 

to first order, since qp is real (see Table I), the imaginary part 
of vpmis negative and proportional to (ka/2)V3. it is assumed here 
that ka is sufficiently large that 
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can be replaced by one. Furthermore when ka sine >> 1, 


(A-18) 


3P 

V, 


szin. 

96 


(-cos) 


% 


H •i'ff 

Ov IT 

e ^ e 


Trka sine 


-JVp0 


+ je 


-j.p(2 


TT-e)j 


• ka 


+ 



2 2/3- 
^p l<a 


From Keller and Levy [15] 


H^^^(ka) 


(A-I9) 


with the soft diffraction coefficient, which is given by Eq. (27) 
after asymptotically approximating the Hankel functions to first order. 
Also 


(A-20) Q H<^>(ka) = (ka), 

Vp N ^.^a Vp 

since H^^^(ka) = 0. 


If we let 


(A-21a,b) .1^^ = 0 - I and - 


and make use of Eqs. (A-18), (A-19) and (A-20), Eq. (A-16) becomes 


(A-22) E -.-^^5!— I 
^ sine p=l 






^p R 


jkR 
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which is an asymptotic approximation for the electric field in the 
shadow region of the H-plane for e < tt. The are the soft EM , 
diffraction coefficients for the sphere given in Table I. 

Next let us use the GTD method described in Section I to 
find E^, By comparing the resulting expression with the asymptotic 
approximation, the launching coefficients and the constant C can be 
found. Consider the field point P in the shadow region of the H-plane; 
according to Keller's generalized Fermat's principle, two surface rays 
emanate from the magnetic current moment at R' = az and traverse paths 
t] and t2 along meri deans to Q-] and Q 2 , where they tangentially shed 
the diffracted rays 1 and 2 toward P, as shown in Fig. 2-A. Thus the 
total far-zone electric field at P is the sum of the fields contributed 
by diffracted rays 1 and 2. It is assumed that ka is sufficiently 
large that the contribution from the surface ray which encircles the 
sphere one or more times before shedding a diffracted ray toward P 
can be neglected. This is consistent with replacing 

• r% 


^ ong in the asymptotic approximation. Since P^ * t' = Pm 
p^ • b = 0, one obtains 

-jks^ “jkS2 

(A-23) = Cp„ G, I- + Cp„ G« I- 

<|) ^m 1 s^ ^m 2 S 2 

wi th 


(A-24a) 

t^ = a(e - 1^) 

(A-24b) 

^2 " ® 

(A-24c) 

S-| ~ S 2 “ R, 

and 


(A-24d) 

- 1 

Jsine 


aij^i , 

) = ai|^2> 

for both and Q^, 
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Upon using Eq. (14), 
CP„ 


(A-25) 


E, = 


m 


y D^le 
^ s/Ii^p=l P P' 


[■ 


-(jk+apa<jji “(jk+ap)aijj2 


+je 


e 
J R 


•jkR 


where the factor j has been inserted to account for the phase jump 
of ir/2 resulting from the caustic on ray path 2. This is the GTD 
expression for E,; comparing it with the asymptotic approximation in 
Eq. (A-22) ♦ 

(A-26) C = - ^ 


(A-27) L= = -. H<^>'(ka)D=. 


(A-28) V = ka - ja^a. 

^ P 

The constant C is chosen to make Eq. (15) consistent with Eq, (17). 
The formula for aS has been obtained previously from o then' canonical 
problems, (see Table I). Since vp ^ ka and ka » 1, h|, 2) (ka) may be 
approximated asymptotically in terms of an Airy function[9] 


,( 2 )' 


-jf, -2/3 


(A-29) (ka) -V - 2e Ai'(-q )1 


' 2 2/3 q„ 


and 

c ~^T9 1/2? 2./3 

(A-30) L® = e ’2(2,k)’/2{^) A1'(-qp) 


1 - fe) 


2/3 


i— 1 

1 n 0^3 

15 V J 


where is the soft EM diffraction coefficient given in Table I. 
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2. E-plane (<t) = |^) Hard EM Boundary 
From Eqs. (A-13) and (A-14) 


(A- 31) 


k^p 


m 


Ir 


-jkR 


7T 

2ka 


i 

n=l 


(-j)' 


2n+l 


H^\^ka) 

n^ 


- Q H^^](ka) 

-JO— 

Q H^^](ka) 

n4 


H^^](ka)f 

n4l J 


‘^Pp^-cose) 

d? 


Eq. (A-31) resembles Eq. (A-15) and we proceed in the same manner to 
a form which can be compared with the GTD solution. Employing the 
Watson transformation and the Cauchy residue theorem as in the 
preceding case 


(A- 32) 


o-JkR , - 2 -Jj 

f-.i me''/TT ve'^e 

^9 J 2— ^ J, '’p >' 

P“' 1 + e P 


Q H^^^ka) 
_E 


'v. r5T — ~ H^^^(ka) P”^ ,(-cose) 

3 ?< u(2)/,_^ v„ ' ' 96 v„ r 


l-QH^^^(ka) > 

9v Vp 


in which Vp are the complex orders for which 

ft H^^Va) = 0. 

(A-33) Vp = ka + (|^)^/^ qp e 

to first order, ^p is real and defined by the above two equations. 
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From Eqs. (A-18) and (A-21), where p in Eq. (A-18) has been 
replaced by its approximation given by (A-33) 


(A-34) 


E 

0 


_ jk ^m y 

sine p=l 


. H(2)(ka)Dj 

'"p P. 


X e ^ + j 




^p R 


■jkR 


which is the asymptotic approximation for the electric field in the 
shadow region of the E-plane for 0 < tt. The is the hard EM dif- 
fraction coefficient given in Table I. P 

The corresponding GTD solution is_deriyed in the ^ame ipanner as 
in the preceding section, except that p^. • t' = 0 and Pm • b = Pm here, 
so that we begin with Eq. (12). The GTD expression for the electric 
field in the shadow region of the E-plane is 


(A- 35) 


CP 


00 



P P 


■ -(jk+a^)aii/i -(jk+a^)aij;2 
e ^ +je ^ 


e-jkR 

R 


9 


Comparing it with Eq. (A-34), 
(A-36) C = - 


(A-37) 


|h _ . /jk^ 

■-p ■ 


1/2 


H^^^ka)o[j, 
> P 


u 

v„ = ka - ja„a. 

P P 

The formula for ap like that of ap, has been obtained previously, see 
Table I. Asymptotically approximating H(2)(ka) in terms of an Airy 
function 


(A-38) 


1/3 


nf^ka) .2 e ^(^) Ai(-q )^1 + 
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and 


u Jt? 1/2. 1/3 
(A- 39) l" = e lM277k) (^) Ai(-qp)x 


1 + 


4 ) 


2/3 


Ts^p ® 




L 

where is the hard EM diffraction coefficient given in Table I. 

The desired launching coefficients are given by Eqs. (A- 30) and 
(A-39), but before we can extend these results to more general convex 
surfaces, it must be determined to what extent they depend on the 
spherical surface employed in the canonical problem. To answer this 
question the radiation from magnetic line sources on a cylinder will 
be treated next. 


3. Cylindrical surfaces 

In the two-dimensional source free region exterior to the closed 
contour L, which is the cross section of the cylinder here, the scalar 
field 


(A-40) ii^Cp) ['l'(p') 1^1 GC^ »p) - GCp* 1^1 (p' )]d£.' 

in which G(F>^' ) is a Green's function, dn' is an infinitesimal distance 
taken along the outward direction to the contour L, and the direction 
of integration is taken in the clockwise sense. The Green's function 
satisfies the partial differential equation 

(A-41 ) [v^ + k^] G(p,p') = - 6{p-p'), 

2 — — 
where is the two-dimensional Laplacian operator and p and p' are 

the position vectors for the field and source points, respectively. 

Eq. (A-40) will be employed to find the electric field radiated by two 

types of axial line sources (see Eqs. (20) and (21) on a circular 

cylinder). 


Let a magnetic line current M = z M 2 be located at p' = a, 

(j)' = 0 on the surface of a perfectly-conducting circular cylinder of 

radius a. On the surface of the cylinder E ^(a,(i)') = M,6(d)')/a. If 

<P ^ 


(A-42) H^(p) = ^(p) 
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and 


(A-43) 1^ = 0 for p' = a. 


r __ 9 H (p') 

(A-44) H^(F) = G^(p,p') -^ - 1 dJl'. 

Since 

8 H 

(A-45) = - jojE = - joie 6((()')/a , 


(A-46) E^(p) = G^(7,;a). 

Next let a magnetic line source M = be located at p' = a, 

(()' = 0 on the surface of a perfectly-conducting surface of radi^^s a. _ 
On the surface of the cylinder E (a,(|)') = - M 6((j)')/a. If E (p) = ij;(p) 
in this case and z <p z 

(A-47) G^(aF,F) = 0, 


(A-48) 4 E^(p') Ij, Gj(? J) in' 




p'=ax 


The Green's functions G|,, Go which satisfy Eq. (A-41), the 
boundary conditions in Eqs. (A-43) and (A-47), and the radiation con- 
dition are 

(A-49) G(p,<|,; a. 0) = - i j 

^ tTFO 


'm 




m 


Q Jn^(ka) 


.( 2 ), 


) 


H;"'(ka)lH;^"^(kp) 


( 2 ), 


m 


in which 


Q = 


±_ 

dp 


for G. 
f h 


and 1 for 


Gg, and 


e = ■> 
m 


1 . 

2 , 


m = 0 
mfO 


cosm(}i 
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Enploying the Watson transformation and computing the residues of the 
integrand again, we transform Eq. (A-49) to a form suitable for an 
asymptotic approximation and ray optical interpretation. 


(A-50) 6(p,4>; a, 0) = 

Q H^^^(ka) 

CO ^ Vp ^ ^ 






v=v_ 


with v„ the zeroes of Q H^^^(ka) = 0, ^We have used 

p V 


sinv. 


1 - e 


-j2-rrv^ 


(A-52a,b) >1^1 = ip 2 = ^ - 4> 

and approximated (1 - e by 1 for ka large. 

Introducing the definition of the diffraction coefficient 


(A-53) 


Q H^^^(ka) 

— ""P 

I^Q H^^^(ka) 

9v Vp 


ik 2 

= - (^) [D„J^ 

^ TT 


op- 


in which Q fo’" Dgp 1 D^p, and the asymptotic approximation 


(A-54) H<2)(kp) .JX 


— -j[kp-(2vfl)|] 

G 3 


the electric fields given by Eqs. (A-46) and (A-48) may be determined 
in the far-zone for ka large. For the circular cylinder D = D . 
Thus . ^ op p 


(A-55) 


-jkM 

■ t I 
s/^kj p=l 





1/2 



] D 


H^^^ka) oh X 

Vp P 

h 
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ro|=j 



and 


(A-56) E = , 

■J 8ttI< j p=l 


-jkM CO 1/2 , 


-jv 4), 

X [e P ^ + 


^ p-jkp 

p qS e_ 


Note that G (p‘ may be obtained from Eq. (A-50) by replacing ^ 
H(2)(ka) in^the^numerator by its derivative with respect to p, kH(2)'(ka), 

P '^P 

Referring to Fig. 2A, the corresponding GTD expressions for E and 
Ez are readily deduced from Eqs. (12), (14), (22) and (24) in 
the manner described in the appendix. Thus 

u h 

“ u u -(jk+a laip, -(jk+a )aij;p _-jkp 
E. = C M_ I L^D^[e P ’+e P 2] S- 


(A-57) 


P = 1 


-P p- 




and 

(A-58) 


pi, 


-(jk+a^)aij;, -(jk+a^)ai|;p -jkp 

'+e ^ . 




Comparing Eqs. (A-57), (A-58) with Eqs. (A-55), (A-56), respectively, 
it is seen that 

i- 

J4 

(A-59) C = - , 

'JS'nk 


(A-60) H^^^ka) D^, 

P ^ 

1 /2 

(A-61) L= = (Jt) H<2)'(ka) (^. 

12 ) It is found. fhat the second order asymptotic approximations of 
H^^^(ka) and (ka) are given by Eqs. (A- 38) and (A-29). which is 
interesting because the vp are different to order (2/ka)4/o in the 
solutions to the hard sphere and hard cylinder problems. Using these 
asymptotic approximations together with the cylinder diffraction coef- 
ficients listed in Table I, one can obtain the hard and soft launching 
coefficients for the cylinder to second order. 
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Comparing Eq. (A-60) with Eq. (A-37) and Eq. (A-61 ) with Eq. 
(A-27) and noting the asymptotic behavior of the Hankel function and 
its derivative just n^ntioned, leads us to conjecture that in the case 
of a general convex surface the launching coefficients are given by 
Eqs. (35) and (36) in the text. 
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APPENDIX II 

THE FIELD IN THE TRANSITION REGION 


In the transition region adjacent to the shadow boundary, the 
representation of the field radiated by the slots described in 
Appendix I is poorly convergent, so a different form is needed. It 
will be seen that a suitable asymptotic solution can be obtained in 
terms of integrals. The first order approximation is usually adequate 
for ka>3, and since the added complexity of the second order approxi- 
mation is difficult to justify, we will derive only the first order 
approximation. Furthermore, since the first order approximations for 
the cylinder and the sphere are the same, except for a factor d^p/d^j 
which is present in the latter case, only the cylinder problem will be 
treated here. 

The Watson transformation is applied to Eq. (A-49) again, but 
this time the resulting expression is left in the form of an integral. 

In the far-zone Hv2J(kp) is approximated asymptotically as in Eq. (A-54) 
so that ^ 


(A-62) 


1 e 


rJkp 


2-J Sirk j yfp z=0 


H^^^(ka)H^^^ (ka)-H^^^ka)H^^^ (ka) 

^ — 

H^'^^ (ka) 


r “jv(<|^1+2ir£) -jv(iJ;/j+ 2 Tr£) 

' e + e 


(A-63) 
8G^(P|Q') 


9p 


T 


K el^'^P ^ 
T— I, 


dv> 


H^^^' (ka)H^^^ka)-H^^^' (ka)H^^^ka) 


2^8Trkj £=0 

-jv(i|;i+2Tr£) -jv(ipp+2TT£)l 

I + e ^ J 


Hj_^^(ka) 


3TT 


where ^ - ir/2, and \p 2 = 2 — ‘5>» Appendix I 
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From the Wronskian relationship for the two types of Hankel functions. 


(A-64) . H^^^'(ka) H^^^ka) - H^^^'(ka) H^^^ka) = 

V V V V TT Ka 


Using Eq. (A-64) in Eqs. (A-62) and (A-63) 


(A-65) 






2i_ 

Trka 




V STikj J 



^00 

H 

- jv ( >|>I +2 ttz) - j v( ijn+2 tt£) 

e + e 

a V 

J-CO 

H^^^'(ka) 

^ V 


and 

(A-66) 


8G5(P/Q') 


9p 


1 'X. 


2jk 1_ el^^P “ 

z=0 




v( l|q+2 iri,) -j 

H^^^(ka) 

V 


v( >1^+2 IT 2,) 


Replacing H(,^^(ka) and (ka) in (A-65) and (A-66) by their Watson 
approximations (forv%ka>>!) gives rise to the Fock function representations 
for and aG /3p'. 


(A-67) G^(P/Q') ^ 


1 


(Sirkj)'’^ £=0 


-Jka(i^i+2TrO ka . e 


T- ^ 
1 ^ 


-jka(tj^2'^2-rrii) 

+ e Y~ 


e‘«2- 

wA(t) 


dx 


e 


-jkp 


and 


(A-68) 

8G3(P/Q') 


3p irka 




(STTkj) 


£=0 


-jka ()|; i+2tt£) . 2/3 

i ^ {^) (-3J7) 


r 

r,W 


d^^ + 


-jka(4,p+2T7i!.) 
+ e 




£_j'^ 

J7 
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wi th 


1/3 

(A-69) H<2>(ka) ■vL (2 ) 

I ■ 2/3 

(A-70) H^2)'(ka) w'(t), 

sj TT 

and 

ka 

(A-71) 5-, = (^) 4/-J ^ 

2 2 * 

Here W2(t) is a Fock-type Airy function related to the Miller-type 
Airy function by 

^ . jf 

(A-72a) W 2 (t) = 24tT e Ai (- re 


alternatively, it may be defined by 


2 

(A-72b) W,(t) = ]= e^^"^ dz 

, tZ-z^/3 

(A-72c) W'(t)=-=L ze di 

with the contour of integration T 2 shown in Fig. 3A. 
Introducing the Fock functions 

(A-73) 

Jtt r-| 2 
and 
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Fig. 3A. Contours of integration for the Airy and Fock functions. 


(A-74) 


O* 1 




r, 


dx, 


(A-75) G^(P/Q') -j ^ 


and 

(A-76) 

where 


-jkaij; 


-JINUl)/. 

9 (^ 1 ) e + 9 (^ 2 ) e 


-jkaij;. 


> S_ 

s/7 


jkp 


9GJP/Q') 2 1/3 

1 


3p' 


C = -jk( 8 TTkj)"^/^, 


"jkai|;-j ^ 

g(?])e + g(?2)e 


-jkaii;. 


, e_ 

r. 


-jkp 


and only the first term £=0 is retained for large ka. When the Fock 
functions are represented by a residue series [16], it is seen that 
they become negligible for c»l, which justifies ignoring the higher 
order terms in the series. 
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Substituting the preceding expressions into Eqs. (A-^6) and 
(^A-48) the far-zone electric fields for the line sources z M and 
(p M, are 

<p 


r -jkaii;, . -aka^p] -jkp 
(A-77) E 'vCMlg(?^)e ' + g(?2^e ^ 

L >/p 


and 


../2 


(A-78) - C 


1/3 


-jkaij/ 


§'(5,)e ■ '+§1(52)6 


■jkai(i. 


e 

Jp 


-jkp 


The above asymptotic representations of E(j) and Ez may be 
interpreted ray optically, as in the preceding appendix. Referring to 
Fig. II-A, the first term in each equation is the field of ray 1, and 
the second term is the field of ray 2. In the transition region 
li^l = ()) -7T/2 is of order (ka)“V3, and for large ka assumed here, this 
angle is small; on the other hand, ii>2 it. Ray 1 traverses only a 
short distance t|i]a on the surface before radiating. The field 
contribution from ray 2 is more accurately represented by the ex- 
pressions of Appendix I, since they are accurate to order 2, and 
this accuracy is important if the ray has traversed a substantial 
distance on the surface between its source and the point of radiation 
Q 2 . However, it is the representation of the field of ray 1 which is 
of interest here, so we will not include the contribution from ray 2 
in the discussion to follow. 


Comparing Eqs. (A-77) and (A-78) with Eqs. (A-57) and (A-58) and 
with Eq. (24) it is seen that 


-jkai|i, 

(A- 79) F = gU^)e ' 


and 


o 1/3 — jkai^ii 

(A-80) G = - j(-j^) g(5-|) e 


These integral representations of F and G may be used in place of the 
poorly convergent series of surface ray modes in the transition region. 
However, the series representations of F and G useful in the shadow 
region can be obtained directly from their integral representations, 
showing that the two solutions join smoothly to first order. This 
will be discussed next. 

In the deep shadow C is large and positive and the Fock integrals 
may be replaced by a residue series[16], so that 
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(A-81) F = 




j— 

Cq e -jka,};, 
^ e 

?pA1(-qp) 


j— 

2 1/3 Jf ^ 

(A-82) G = - j(y e ^ I ^ ® 


1 


p=l 


■ATT^ 


which equal 


(A-83) F = 


u u -(jk+a )ai|i, 

I e P ' 

A- n n 


P P 


p. c -(jk+a^)ai|j, 

(A-84) G = I l"d^ e P 


■P P 


when the first order approximations for the diffraction coefficients 
and attenuation constants, see Eqs. (26) through (28), and the first 
order approximations for the launching coefficients 


, -JT2 1/2 2 

(A-85) L=p=e 


(-q„)D^^ 

' ^p' op 


. 1/2 2 1/3 _ u 

(A-86) (|j) Ai(.,p)Djp 


are employed. The two equations above follow from Eqs. (A-27) and 
(A-39). 

It is also important that the integral representations of the 
field in the transition region join smoothly with the geometrical 
optics field, given by Eqs. (18) and (24), in the illuminated region. 
Deep in the illuminated region ^ is large and negative with the result 
that the Fock functions may be approximated asymptotically [16]. 


68 



(A-87) g(d ''' 2 e 


J5^/3 


(A-88) g(5) - 2j5ej^ . 

Substituting the above expressions into Eqs. (A-79) and (A-80) 

-jkasintj;, 

(A-89) F -V 2 e ' 

-jkasini(/, 

(A-90) G 'X. -2 e ' 

3 

in which we have employed {p^ - <p^ /6 sin^j;-!. 

Equation (A-89) is recognized as the geometrical optics form of 
F, when one notes that the exponential term gives the phase delay 
between the point of radiation Q-| and the source located at Q' , see 
Fig. 4A. As pointed out in Chapter III, in the illuminated part 
of the transition region, one visualizes the surface ray as traveling 
from the source at Q' to the point of radiation at Qi where it sheds 
a ray tangentially toward P as indicated in Figs. 7 and 4A. The 
representation of F given in Eq. (A-79) may be used in both the 
illuminated and transition regions. 



Fig. 4A. An equivalent ray system for calculating 
the radiation in the illuminated region. 
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On the other hand, the asymptotic form of G in Eq. (A-90) does 
not join smoothly with the geometrical optics solution, 

-jkasini/», 

(A-91 ) G - 2 sini|)^e , 


unless one modifies Eq. (A-80) in the following way 


(A-92) 
wi th 


2 1/3^ jkail^. 


f 


(A-93) fU>,) = 


simf/-! 
1 


- -2 < 'J'l 1 0 


^ 0 


Equation (A-92) has been justified through its application to 
numerous examples; it should be regarded as an algorithm for 
calculating G in the transition region so that it joins smoothly 
with the geometrical optics field, which is valid deep in the 
illuminated region. 
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APPENDIX III 

DERIVATION OF THE FIELDS AT AN AXIAL CAUSTIC 
OF THE DIFFRACTED RAYS 


Consider a perfectly-conducting convex surface of revolution 
with an infinitesimal slot positioned on its axis of revolution at 
Q' as shown in Fig. 5A. There is a caustic of the diffracted rays 
on the axis of revolution in the shadow region and the ordinary GTD 
solution breaks down here. This appendix describes how the GTD 
solution may be supplemented by an integral representation of the 
field to obtain a high frequency approximation at the axial caustic. 
However, it will be seen that the GTD is used to deduce the 
equivalent sources of the integral representation. The field in the 
shadow region (excluding the negative z axis) is obtained by a super- 
position of the fields associated with the rays shed from only two 
points Q] and Q 2 on the surface of revolution, as shown in Fig. A2. 

In the vicinity of the caustic (z-axis), the rays arriving at P are 
shed from a continuum of points constituting a ring of radius 'a' on 
the surface as shown in Fig. 5A; the radial distance 'a' is the maximum 
perpendicular distance from the z-axis to the surface. 




(CAUSTIC OF DIFFRACTED RAYS) 

Fig. 5A. Ray geometry for calculating the field at an axial caustic. 
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The equivalent source for the infinitesimal slot at Q' is a 
magnetic current moment of strength pmX at Q' which is oriented along 
the slot axis as shown in Fig. 5-A. This equivalent source produces an 
electric field at the point, T located at p=a, 0=<j>' , z=-Zq. This point 
is in the near zone of the surface of revolution, but it is several 
wavelengths from Qq, the point at which the surface ray sheds toward T. 

The GTD result for the electric field at T is 


_ ^ ^ I dip 

(A-94) E(T) = C tnsinVF„ - ^ • 


where 

(A-95) 


I L^(Q')D^(Qo)e ^ 

p=l ^ ^ 


f 0 f 0 h 

-jk dt' - a"(t')df 
Jn* Jn* P 


-jk 


(A-96) = 






dt' 

Q* o ^Q* 


a„(t')df 


C = - jk/4TT 


and the quantity dn, and the caustic distance, p are defined in 
Chapter II. Note that Fq and Gp ar^ the same as the F and G 
functions without the factor^ dl^o/d'i'. For the ray shed in the 
caustic direction, p->«>, so that Eq. (A-94) becomes 


1 

(A-97) E(T) [C Pj^{nsin<t)'F^ + 6coS()>'G^} -!=r ] •— 

U Uo 

The integral representations for the electric fields at T radiated 
by electric and magnetic ring currents of strengths I (cf)"), and 
M ((;>"), respectively are 


(A-98) 


Ee(T) 


jka 


4ir 


2tt 




[Zg(rxrxl^(<j>")4)")] 


,-jkr 


d<j>". 


72 



and 


(A-99) 


"m 


(T) ^ 


jka 

4tt 


•2ir 


[rxM^ (<!)")<})"] 


,-jkr 


d<|)" 


where 



a is the radius of the ring currents. 


r 



+2a^[l -cos (<}.'-<()")] 


5 


and the subscripts e and m refer to the electric field contributions 
from the electrig type and magnetic type ring currents, respectively. 
The unit vector r is directed from (|)" on the ring to T. If one 
evaluates the above integrals by the method of stationary phase, one 
finds that there is a stationary phase point at Qq. The contributions 
to the field at T which arise from this stationary point are 


J’t 

(A-lOO) 


and 


(A-101) 

~o 

A comparison of Eg. (A-97) with Eq. (A-lOO) and Eq. (A-101) indicates 
that the GTD field at T may be considered to have the equivalent 
electric and magnetic ring currents whose strengths Io(<|)') and M ((j)') 
respectively are given by 




4 


(A-102) 


-jfo 1/2 

C Pm ^0 ^0 ® = 


73 



and 


(A-103) 


Mq(<|.') = - 


C sin<|,' 




1/2 


The radiation field in the vicinity of the caustic is therefore the 
field radiated by these equivalent ring currents. This technique of 
using equivalent currents based on the GTD to describe the fields in 
the vicinity of the axial caustic of the diffracted rays has been used 
with success elsewhere in problems involving curved edge diffraction [17,1 8]; 
however, unlike the curved edge, the rim indicated in Fig. 5A does not 
stay fixed as the observation point moves slightly away from the caustic, 
and hence the caustic solution presented here is not expected to be 
accurate at aspects not close to the axial gaustic. The expression 
for the far-zone electric field due to the (})' directed electric and 
magnetic ring current sources of strengths Io(^') and MgCcj)') given by 
Eqs. (A-102) and (A-103), respectively is expressed in terms of the 
radiation integral as 


(A-104) 


jka e 


E0>) 


-jks 


'2tt 


[Z (fx?xl (♦: )}■ )+(raM„{(,' )}' )]eJI<asinecos( 4 ,-V _ 


This expression for E(P) in Eq. (A-104) is a general result, and the 
integrals involved can be evaluated in closed form; however only the 
integrations for the patterns in the X-Z and the Y-Z planes are needed, - 
see the discussion connected with Eq. (60) of the text. 

X-Z plane (((>=0) 

(A-105) 

^ Pm [6Q{jQ(kasine)-J2(kasine)}+F^cose{Jjj(kasine) + 

- j ks 

+ J2(kasine)}] ^ — 
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Y-Z plane ((f)=^) 

(A-106) 

1/2 

Ej^(P) ^ [GQC0s6{J^(kasine)+J2(kasine) }+Fp{jQ(kasin9) 


" j ks 

- J2(kasine)}] — 

in which Fg, Gq are given by Eqs. (A-95), (A-96). The JqCx) and J 2 (x) 
in Eq. (A-105) and Eq. (A-106) are cylindrical Bessel functions of the 
first kind with argument x(=kasine); the subscripts indicate the orders 
of these functions. 

As mentioned before, the above results are valid only when P is 
in the vicinity of the axial caustic. This means that although the 
caustic correction is needed for the angular range tt- (ka)"V3 < e _< tt, 
the above formulas may be satisfactory only for a smaller angular 
range adjacent to tt. This is particularly true when ka is small. 

Work is currently under way to improve these formulas so that they 
are adequate throughout the entire range tt- (ka)"V3 < e tt. 
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